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Abstract 

It is pointed out that string-loop modifications of the low-energy matter 

couplings of the dilaton may provide a mechanism for fixing the vacuum 

Q>^ ■ expectation value of a massless dilaton in a way which is naturally compatible 

^«0 , with existing experimental data. Under a certain assumption of universality 

of the dilaton coupling functions , the cosmological evolution of the graviton- 

CP \ dilaton-matter system is shown to drive the dilaton towards values where it 

decouples from matter ( "Least Coupling Principle" ) . Quantitative estimates 

are given of the residual strength, at the present cosmological epoch, of the 

coupling to matter of the dilaton. The existence of a weakly coupled massless 

Qh! dilaton entails a large spectrum of small, but non-zero, observable deviations 

(— I I from general relativity. In particular, our results provide a new motivation 






X 



for trying to improve by several orders of magnitude the various experimental 
tests of Einstein's Equivalence Principle (universality of free fall, constancy 
of the constants,. . . ). 

I. INTRODUCTION 

At present we know only one theory which treats gravity in a way consistent with quan- 
tum mechanics: string theory. In the low energy limit (low in comparison with the Planck 
mass) string theory gives back classical general relativity, with, however, an important dif- 
ference. All versions of string theory predict the existence of a (four-dimensional) scalar 
partner of the tensor Einstein graviton: the dilaton. It may happen that this scalar field 
acquires a mass due to some yet unknown dynamical mechanism. This is the generally 
adopted view, and if so there will be no observable macroscopic difference between string 
gravity and Einstein gravity. In this paper we will discuss another possibility: that the dila- 
ton remains massless. This immediately leads to the dramatic conclusion that all coupling 
constants and masses of elementary particles, being dependent on the dilaton scalar field, 
should be, generally speaking, space and time dependent, and influenced by local circum- 
stances. This conclusion is of course not new and it was precisely the reason for discarding 



the possibility that we are going to discuss. Indeed, it has been stated that the existence of 
a massless dilaton contributing to macroscopic couphngs would, at once, entail the following 
observable consequences: (i) Jordan-Fierz-Brans-Dicke-type |I|] deviations from Einstein's 
theory in relativistic [0{Gm/c^r)] gravitational effects [§]; (ii) cosmological variation of the 
fine structure constant, and of the other gauge coupling constants 0, and (iii) violation of 
the (weak) equivalence principle Q] . As the strength of the coupling of the dilaton to matter 
is expected to be comparable to that of the (spin 2) graviton, and even larger than it in the 
case of hadrons ^, the above observable consequences seem to be in violent conflict with 
experiment. Indeed, present experimental data give upper limits of order: (i) 10~^ for a pos- 
sible fractional admixture of a scalar component to the relativistic gravitational interaction 
IP, (ii) 10"^^yr~^ for the fractional variation with time of the fine-structure constant^ P|, 
(iii) 10~^^ — 10~^^ on the universality of free fall (weak equivalence principle) |f^-p!0| 0. [See 



rT| , p!2[| for reviews of the comparison between gravitational theories and experiments]. 

In this paper, we point out that non-perturbative string loop effects (associated with 
worldsheets of arbitrary genus in intermediate string states) can naturally reconcile the 
existence of a massless dilaton with existing experimental data if they exhibit the same kind 
of universality as the tree level dilaton couplings. By studying the cosmological evolution of 
general graviton-dilaton-matter systems we show that the dilaton is cosmologically attracted 
toward values where it decouples from matter, a situation which we call the "Least Coupling 
Principle" . Roughly speaking, the origin of the attraction is the following. Masses of different 
particles depend on the dilaton, while the source for the dilaton is the gradient of these 
masses. It is therefore not surprising to have a fixed point where the gradient of the masses 
is zero. [With some important differences discussed below, this mechanism is similar to the 
generic attractor mechanism of metrically-coupled tensor-scalar theories discussed in Refs. 
|r3|]. This cosmological attraction is so efficient that the presently existing experimental 



limits do not place any significant constraints on the physical existence of a massless dilaton. 
Most importantly, we give quantitative estimates for the level of residual deviation from 
Einstein's theory expected at the present cosmological epoch, notably for the violation of 
the equivalence principle. 

II. THE GRAVITON-DILATON-MATTER SYSTEM 

At the tree level in the string loop expansion (spherical topology for intermediate world- 
sheets) the effective action describing the massless modes (here considered directly in four 
dimensions) has the general form |]l^|-p!6[ 



^tree = / d^X^QC-^'' {{a'Y^R + 4 □$ - 4(V<I')^ 



^Note that, within the QCD framework, it does not make sense to speak of the variation of any 
strong-interaction coupHng constant (the hadron mass-scale adjusting itself such that Ostrong — 1). 

^The most recent analysis of Lunar Laser Ranging data |jlO| finds that the fractional diff'erence 
in gravitational acceleration toward the Sun between the (silica-dominated) Moon and the (iron- 
dominated) Earth is (-2.7 ±6.2) x 10^^^ 



-'^F^^F^^'' - ^D^ + 



+ j:0[ia'd')-]}. (2.1) 



n>l 



Here, c/f^u (often denoted G^i,) denotes the metric appearing in the cr-model formulation 
of string theory and is used for defining all the covariant constructs entering Eq. (1) 
[V, (-F)^, D,- ■ ■]; $ denotes the dilaton; a summation over the various possible gauge fields 
[F;^ = d^Al - d^A"^ + P^^A^Al] and fermions \D = Ti^t. + Ap'')] is understood; the 
ellipsis stand in particular for the ill-understood remaining scalar sector of the theory [Higgs 
fields and their Yukawa couplings, and possibly other gauge-neutral scalar (moduli), or 
pseudo-scalar (axion,. . . ), fields]; and the last term symbolically denotes the infinite series 
of higher- derivative terms representing the low-energy effects of all the massive string modes 
on which one has to integrate to get the effective action for the massless modes. 

The remarkable feature that, when formulating the action in terms of the "string frame" 
metric Qfj,^, the dilaton couples, at the string tree level, in a universal, multiplicative manner 
to all the other fields derives from the fact that Qs = exp(<l>) plays the role of the string 
coupling constant. In the cr-model formulation, this is easily seen to follow from applying the 
Gauss-Bonnet theorem [{Atc)~^ J d'^^^/hR^'^\h) = x = 2(1 — n); n= number of handles] to 
the Fradkin-Tseytlin a-model dilaton term, Sdn = {4:7r)~^ J (P^y/h^{X)R^'^\ In a constant 
(or slowly varying) dilaton background, the genus-n string-loop contribution to any string 
transition amplitude contains the factor exp(— ^du) = exp(2(n — 1)$) = (/g*^"^^-*. Therefore, 
when taking into account the full string loop expansion, the effective action for the massless 
modes will take the general form 






-5i.($)^F2-5^($)jD^+...| . (2.2) 

At this stage of development of string theory, one does not know how to control the structure 
of the various dilaton coupling functions -Bj($) (i = g,^,F,ip, ...) beyond the fact that in 
the limit $ -^ — oo {gs — > 0) they should admit an expansion in powers of g"^ = exp(+2$) 
of the form, 

5,($) = e-2* + 4^^ + cPe^* + c?e^* + • • • (2.3) 

[Note that we have in mind the low-energy regime , with broken supersymmetry, for which 
there are no a priori obstacles to having couplings of the type ( |2.3| ) with c^*^ 7^ 0.] 

Concerning the low-energy effects of all the massive string modes, we shall assume for 
simplicity that, like at tree level, Eq. (PTT|), they are equivalent to introducing a cut-off at a 
^-independent string mass scale A^ ~ (a')~^/^, when measuring distances by means of the 
string-frame metric g^j,,^. 

It is convenient to transform the action ( p.2| ) by introducing several ^-dependent rescal- 
ings. One can put both the gravity and the fermion sectors into a standard form by: (i) 
introducing the "Einstein metric" , 



g^, = C Bg{^)g,,, (2.4) 

(with some numerical constant CQ), (ii) replacing the original dilaton field $ by the variable^ 

1/2 

(2.5) 



^-h'\( 



^^^iV + 2^ + 2^ 

Bg) Bg Bg 



(where a prime denote d/d^), and (iii) rescaling the Dirac fields 

^ ^ C-'/%'/^B]!'^ . (2.6) 

The transformed action can be decomposed into a gravity sector (^f^j,, ip) and a matter one 

S[g,ip,'ip,A,---] = Sg^^ + Sm , (2.7a) 

Sm = ld''x^l-^D^-jBF{ip)F^ + --\ . (2.7c) 

Here, g = AnG = \Ca' {G denoting a bare gravitational coupling constant), Bp{ip) = 
Bp[^{ip)] and the ellipsis stand for the (more complicated) Higgs sector. One should note 
that the string cut-off mass scale acquires a dependence upon the dilaton in Einstein units: 

A,{^) ^ G-^'^B-"\^)K . (2.8) 

Essential to the following will be the dilaton dependence of the matter Lagrangian. One 
does not know at present how to relate string models to the observed particle spectrum. The 
basic clue that we shall follow is the dilaton dependence of the gauge coupling constants: 
g'"^ = kBp{ip) from (|2.7c| ). To connect the (bare) effective action (|2.7] ) (integrated over 
the massive string modes) to the low-energy world, one still needs to take into account the 
quantum effects of the light modes between the string scale As(v9) and some observational 
scale. In the case of an asymptotically free theory the ratio of the IR confinment mass scale 
Aconf to the cut-off scale, is, at the one-loop level, exponentially related to the inverse of the 
gauge coupling constant appearing in the bare action: 

Aconf ~ As exp(-87r2r 1(7-2) = G-^'^Bg^/^^) exphSTr^r^A; 5^(<^)]A, , (2.9) 

where the one-loop coefficient h depends upon the considered gauge field as well as the matter 
content. The mass of hadrons is, for the most part, generated by QCD-effects and is simply 



^We shall choose C such that the string units and the Einstein units coincide at the present 
cosmological epoch: Ci?g(<I>o) = 1- 

^When Bg = B^ the quantity under the square root in Eq.(P^) is positive definite. 



proportional to Aqcd ( with some pure number as proportionality constant). The dilaton 
dependence of the QCD part of the mass of hadrons is therefore given by (|2.9| ) with 6 = 63 
and Bp = B3 being the appropriate QCD quantities. However, the lepton masses, and the 
small quark contributions to the mass of hadrons, are not related to Aqcd ( at least in any 
known way). Their dilaton dependence is defined by specific mechanisms of spontaneous 
symmetry breakdown (and compactification) which depend on particular string models and 
are not well established at present. Let us note that in technicolor-type models, as well as in 
no-scale supergravity ones, all the particle mass scales are related to the fundamental cut-off 
scale by formulas of the type (|2.y|). 

As a minimal ansatz, we can assume that the mass (in Einstein units) of any type of 
particle, labelled A, depends in a non-trivial way on the VEV of the dilaton through some 
of the functions Bi appearing in (p.2|): 



rriAi^) = rriAlBgi^), Bf{^), ■ ■ ■] . (2.10) 

The essential new feature allowed by nonperturbative string-loop effects (i.e. arbitrary func- 
tions i?j($), Eq. ( |2.3| )) is the possibility for the function 771^(^9) to admit a minimum for 
some finite value of ip. Assuming this, we shall see below that the cosmological evolution 
naturally attracts (f to such a minimum. However, if the various coupling functions Bi{(f) 
differ from each other the minima of mA{<f) depend , in general, on the type of particle 
considered. It will be seen below that this weakens the attraction effect of the cosmological 
expansion, and, more importantly, leaves room for violations of the equivalence principle at 
a probably unacceptable level (see the footnote following Eq.( |6.13| )). This suggests to con- 
centrate on the case where string-loop effects preserve the universal multiplicative coupling 
present at tree-level, Eq.(^.lD, i.e. the case where all the dilaton coupling functions coincide: 
Bi{ip) = B{ip) for i = F,g,- ■ ■. In this "universal B{ip)" case, the extrema of the function 
rriAif) — ''^a[B{(p)] will (generically) coincide with the extrema of the function B{ip). As 
discussed below, this assumption leads very naturally (without fine-tuning, or the need to 
inject small parameters) to a situation where the present deviations from general relativity 
are so small as to have escaped detection. When we shall need in the following to estimate 
quantitatively the dependence of particle masses on (p, we shall assume that the mass of any 
particle A is of the form suggested by Eq.(f^.9D: 

rriAiip) = fiAB-^'\^) exp[-87rVAS(y;)]A, , (2.11) 

with Ha and va pure numbers of order unity. We believe that our main quahtative conclu- 
sions do not depend strongly on the specific form of the assumption ( |2.11| ). 

For the quantitative estimates below we need to choose some specific value of the string 
unification scale A^ oc a'"-*^/^. The theoretical value A^ = e'^^~'^^/^3~^/^5'sAfpianck/47r ~ (7^ x 
5.27 X lO^^GeV has been suggested fl^. Here Qs denotes the common (modulo possible 
Kac-Moody level factors of order unity) value of the gauge coupling constants at the string 
scale. To fix ideas, we shall take A, = 3 x lO^^GeV. 



III. CLASSICAL COSMOLOGY WITH A DILATON 

The gravitational field equations derived from Eqs. (|2.71 ) read 



R^^ = 2d^ifd^ip + 2q\T^^--Tg^^\ , (3.1a) 

Uip = -qa , (3.1b) 

where the source terms, defined by T^'^ = 2g~^/'^5Sm/8gi_iv, cr = g~^^'^5Sm/S(p, are related by 
the energy balance equation: V uT^^ = aV^ip. 

In the case of a Friedmann cosmological model, ds'^ = —dt^ + a?{t)d(? with dl"^ = 
(1 - Kr^y^dr^ + r^{d9^ + siii^Odip^), K = 0, +1 or -1, the field equations give (T/;' = 
diag(— p, P,P,P); H = a/ a, the overdot denoting d/dt) 

-3- = g(p + 3P) + 202 , (3.2a) 

a 



K 



3H' + 3— = 2qp + if' , (3.2b) 



(p + 3Hip = qa. (3.2c) 



In the following, we concentrate on the spatially fiat case {K = 0). Following Ref. |]T3 
we can combine Eqs. (|3.2|) to write a simple equation for the cosmological evolution of the 
dilaton with respect to the logarithm of the cosmological scale factor: p = ln(a) + const, 
(not to be confused with the pressure P). Denoting d/dp by a prime, one gets {K = 0) 

^" + (1 - A)^' = - , (3.3) 



3-v?'2" ^ '^ p 

where A = P/p. 

Except during phase transitions, the material content of the universe can be classically 
described as a superposition of several (weakly interacting) gases labelled by A, i.e. by an 
action of the form 

Sm[g,^p,XA] = -Y. mA[ip{xA)][-g^^u{x\)dx'Xdx''AY^'^ (3.4) 

A -^ 

(the massless particles being obtained by taking the limit rriA — * with ttiaUa = ruAdx^/dsA 



fixed). In Eq. (|3.4|) the summation over A includes a sum over the statistical distribution of 



the A-type particles. The gravitational source terms corresponding to Eq. ( |3.4| ) read 

r^"'(x) = Yl I dsA'mA['p{xA)\u^A'^''A^^'^\x - xa) , (3.5a) 

.ln(x) A •' 

1 /■ 

^=Y / dsAaAVpi.^A)]mA{'p{,XA)]5^^\x - xa) 



J9[x 
aix^ 



/g{x) A 
J2aAMx)]TA{x) , (3.5b) 



where 



_ din mAJip) , . 

«a(^) = Q- (3.6) 

measures the strength of the couphng of the dilaton to the A-type particles. In the second 
Eq. ( p.5b| ) Ta = —pA + 3Pa denotes the trace of the A-type contribution to the total 



T'^'^ = J^aTa" ■ It is easy to see that when the different A-gases are non interacting their 
corresponding sources satisfy the separate energy balance equations: VuTa" = cfaV^^p = 
aATAV^ip. 

In the string context, it is natural to assume that the string scale A^ subsumes both 
what is usually meant by "Planck scale" and "GUT scale" , leaving essentially no room for 
a quasi-classical inflationary era. We leave to future work a discussion of primordial stringy 
cosmology, and content ourselves by describing the evolution of dilatonic cosmologies through 
a radiation-dominated era, followed by a matter-dominated one. 

IV. EVOLUTION OF THE DILATON DURING THE RADIATION-DOMINATED 

ERA. 

During a radiation-dominated era (universe dominated by ultra-relativistic gases) the 
gravitational source terms are approximately given by 

p^3P^^.(T)^T\ (4.1a) 

a ^ , (4.1b) 

where g*{T) = J2Bose9AiTA/T)^ + {7/8)J2Fermi9AiTA/T)'^ is the effective number of rela- 
tivistic degrees of freedom in the cosmic soup at temperature T. [ The sum defining g*(T) is 
taken only over particles with mass ttia <^ T; because of possible previous decouplings the 
corresponding relativistic gases may not all have the temperature T, e.g. T^, = (4/ll)^'^T^ 
below 1 MeV ]. Eq. ( [4.1b| ) suggests that the dilaton does not evolve during the radiation 



era. More precisely, Eq. ( p.3|) with A ~ 1/3 shows that (f{p) behaves as a particle, with 
velocity-dependent mass, submitted to a constant friction. In a few p-time units, if{p) will 



exponentially come to rest, [see Ref. |]T3| for the exact solution of the damped evolution of 
<f{p) when cr/p is negligible]. However, something interesting happens each time the uni- 
verse cools down to a temperature T ~ rriA defining the threshold for the participation of 
the species A to the relativistic soup. When T ~ rriA, the term on the right-hand side of 
the p-time evolution of (p is well approximated by 



^r = ~:^ —rfA'r±i^A)aA{^) , (4.2) 

Ptot T^ g*[-i- ) 



where za = rnA/T and 



2 r°° , [x^ - z^^'^l'^ 



r±{z)^z'l dx ' ^^^'^ , (4.3) 



where the upper (lower) sign corresponds to A being a fermion (boson). In the approximation 
(justified by the results to be discussed) where the dilaton contributions to the Einstein 



equations ( p.2a]) , ( p.2b| ) are negligible one has T oc a ^, and therefore p = \nzA (with 



an adapted choice of origin for p). Then, as a function of p, the (everywhere positive) 
function t± is proportional to exp(+2p) when p — >■ — cxo, rises up to a maximum ~ 1.16 
when za — 0.87, and falls quickly to zero as exp[|p — exp(p)] when p — > +00. [ This 
maximum occurs because for T ^ rriA the (ultra-relativistic) particles do not contribute to 
a, while for T ^ rriA there are exponentially few particles ]. Remembering the definition 
( ^.6[ ) of aA{^), it is easy to see that if the function mA{^) has a minimum, say ip^, and if 
the initial value of (p, say ip^ = ip{p = — cxo) = ip(T :$> m^), is sufficiently near ip^, Eq. ( p.3|) 
will describe a damped, transient nonlinear attraction of ip{p) around (p^. [Note that, from 
the above discussion, the initial velocity is zero to an exponential accuracy ~ exp{pA' —Pa)]- 
The existence of such an attraction mechanism by mass thresholds during the radiation era 
was noticed in Ref. [^ in a related context (generalized Jordan-Fierz-Brans-Dicke theories 
characterized by a universal, A-independent coupling function ^^(v^) = c^if) = da{ip) / dip) . 
In the context of Ref. |TB[ , it seemed natural to assume that the curvature of the function 
\Y\.m{(p) = a{(p) + In mo near its minimum was of order unity. This rendered the presently 
discussed attraction mechanism very ineffective. An important new feature of the present, 
dilatonic, context is that the curvature of In 777.^(9?) near its minimum is expected to be large 
compared to one. This follows from the expected exponential dependence on Bp{ip) of the 
mass scales of the low-energy particle spectrum, Eq. ( p.9| ). To fix ideas and be able to make 
some quantitative estimates, we shall take the form (|2.11|) with /x^ = 1. This yields 



oia{^) 



lu^ + i 

niA 2 



^MM = +inlL^M:M, (4.4) 

dip niA dip 



where A', = e^/'^Ks ^ 5 x lO^^GeV. 

We see that a minimum (pm of ?77^((y9) corresponds to a maximum of B{ip) (or a minimum 
of B~^{ip)). Let us denote by k the curvature of the function \n.B~^{ip) near its minimum 
ipm- In the parabolic approximation 

\nB'\p^) ~ \nB'\p^,n) + \^^{^ - ^mf , (4.5) 

one gets 

aA{'P>) = (3a{'P> - ^m) , (4.6a) 

(3a = K ln(A;/mA) = k[A0.75 - ln(mA/lGeV)] . (4.6b) 



Inserting Eq. ( [4.6a| ) into Eq. ( [1.2| ) and then into Eq. ( p.3| ) (written in the approximation 



A ~ 1/3) yields 

ip"{p) + ip'ip) = sa{p) b(p) - ^m] , (4.7) 

with 

SA{p) = -^jA^T^{en. (4.8) 

^T^ a*(T) 



Within a good approximation one can replace the temperature-dependent quantity gA/g*(T) 
by its initial value, say /^ = QaI Q^^i in which gf™ = g^{T ^ m^) contains the contribution 
7gA/8 (or qa) if A is a fermion (or boson). Eq. ( ^.7| ) describes a damped motion submitted 
to a transient harmonic force tending to attract if toward ipm- The final outcome of this 
motion is to leave (when p = +00) ip nearer to ipm than it was when it started at rest 
at p = — cxD. We define the attracting factor of the A-th mass threshold as m-|-(6yi) = 
{ip{+oo) — ip^)/{ip{— 00) — v^m)' where the suffix ± in the left-hand side corresponds to the 
fermion/boson case and where Ba = PaIa — PaQa/qT- There are two quite different regimes 
in this mass-threshold attraction mechanism: when 6^ -C 1 (6^ < 0.5 sufficing), f{p) moves 
monotonically toward ipm by a small amount given by integrating over p the force term on 
the right-hand side of ( [4.7|) evaluated at the original position of ip ("kick" approximation). 
The result is (see Ref. [ p!3[] ) 



m^ibA) = I - -bAi +0(6^), (4.9) 




where the upper (lower) coefficient corresponds to the fermion (boson) case, respectively. In 
this first case the attracting power of the A-threshold is rather weak (hence the conclusion 



of Ref. |T3[ that the total radiation era attraction is rather ineffective in the case of usual 
tensor-scalar theories with Pa = 0(1) and Sa/a ~ ln(100/10) ~ 2.3). By contrast, in the 
present, dilatonic context one expects k ~ 1, /3^ ~ 40 and therefore 6^ 3> 1 for many mass 
thresholds (the most efficient mass thresholds being the latest in the radiation era which 
tend to have the largest /™'s: notably the e~^e~ threshold with /™ = 4/10.75 ~ 0.372). 
When 6a ^ 1 (6a > 2 sufficing in practice) one can analytically solve Eq. ( [4.7[ ) by a WKB- 
type approach. [With some subtleties compared to the usual WKB approximation as the 
matching between the damped and oscillating regions must be done via Bessel functions 
instead of the usual Airy ones]. Qualitatively the motion of v^(p) begins by a slow roll 
toward ipm, continues by WKB oscillations around ipm, and terminates as a damped inertial 
motion. The final analytical results for the attraction factor reads (6a ^ 1) 

m±(6A) = (C±6a)"^/^cos^^ , (4.10a) 

with C+ = 15/8, C_ = 15/4 and 

Ot = r^[-SA{p)Y''dp ~- = bfl^ - ^ , (4.10b) 



with /+ ~ 1.2743, /_ ~ 1.4029. Note that when 6a — >■ 00, |m±(6A)| tends to zero as 0(6^ ). 
Fig. 1 represents the two functions 'm±{b), obtained by numerically integrating Eq. ( [4.7|) . 

One must take into consideration the fact that mass thresholds can occur only for parti- 
cles whose masses are smaller than the critical temperature of the phase transition through 
which they acquired a mass (e.g. the pions are the only hadrons to take into account). 
The Higgs threshold is to be considered as part of the electroweak phase transition, and 
the strange quark threshold overlaps with the quark-hadron phase transition. This leaves 
nine, clearly present, mass thresholds associated (in decreasing temperature scale) with the 
top quark (/f = 12/106.75, A ~ 35.74k), the Z^ {f^ = 3/95.25, pz ^ 36.24k), the W^ 
(f^ = 6/92.25, pw ^ 36.37k), the bottom quark (/^° = 12/86.25, pb ~ 39.14k), the tau 



(/^ = 4/75.75, Pr ^ 40.17k), the charmed quark {f^ = 12/72.25, pc ^ 40.35fi;), the pions 
(/;" = 3/17.25, p^ - 42.74fi;), the muon {fj^ = 4/14.25, P^ ~ 43.00k), and the electron 
(/i" = 4/10.75, pe ^ 48.33k). The quoted values of /jf and Pa show that 6^ = Pa/a is 
typically a few times k (with extreme values 1.14k and 17.98k for the Z and e respectively). 
A look at Fig. 1 shows immediately that if k is of order unity, each mass threshold will a be 
rather efficient attractor. The compound effect of all those attr actors is discussed below. 

Besides mass thresholds, phase transitions provide another possible attractor mechanism 
for the dilaton during the radiation-dominated era. During a phase transition the vacuum 
energy density V changes from some positive value, say V^^ = 5'vac(7'"^/30)T^'^, when T > Tc to 
a comparatively negligible value when T < T^. For instance, in the case of the QCD (quark- 
hadron) phase transition one has Tc ~ 200 MeV and g^^c = 34/3 (in a simple model |l^ 
describing the unconfined phase as a relativistic gas of gluons and u and d quarks — besides 
7, e, u and /i — and the confined phase as a relativistic gas of pions). Besides its dependence 
on the temperature the vacuum energy density is also a function of the dilaton. Therefore 
the vacuum term in the matter action, S'vac = —^/9V{ip,T), will generate a corresponding 
source term cXvac = —dV/dip in the right-hand sides of the dilaton evolution equations ( p. lb] ), 
( p.2c| ) or ( p.3| ). In the simple model of the QCD phase transition just described, one can 



estimate the source term cXvac by assuming that the dilaton dependence of V is essentially 
contained in the yj-dependence of the critical temperature Tc. In turn, the latter dependence 
is obtained from Tc ~ ^qcd with Aqcd(v^) given by Eq. ( p.9|) with the appropriate one- 
loop coefficient. This shows that if will be attracted toward a minimum of Aqcoif)- More 
precisely, if we assume, to fix ideas, that Bg{ip) = Bp^ip) in Eq. (|2.9| ) and that ip is near 
the maximum (pm of B{(p), one gets, in the parabolic approximation ( |4.5| ) (setting as above 
A = Ptot/Ptot and p = \n{Tc/T) ) 

^"{p) + ^[1 - X{p)]v'ip) = S.ac(p)b(p) - fm] , (4.11) 

where Svac(p) - -6/9vac/vacexp(4p) when p -^ -oo, with P^^c = Kln(A(,/Tc) and /vac = 
gvac/g*(T > Tc). After the phase transition, when p — > +oo, one expects Svac(p) to fall 
quickly to zero as exp(— aT^/T) = exp(— aexp(p)) with a of order unity. In the limit where 
6vac = /^vac/vac IS large enough to make ip oscillate around ipm, one can solve Eq. ( ^.11| ) by 
a WKB-type approach. The final result for the attraction factor due to a phase transition, 
p(&vac, ■ ■ ■) = ifi+'X)) - ipm)/i<p{-oo) - (fm), reads 

J9(&vac, • ■ ■) = 23/27r-ir(5/4)ai/2(65^^J-i/8exp(_/) cos^ , (4.12) 

where / = (1/4) J^^[l — 3\{p)]dp, and where the angle 9 depends on the two functions \{p) 
and Svac(p)- [In the approximation X(p) = 1/3, one finds 9 = J^^[—s^g,cip)Y^'^dp — n/8]. 
In the case of the QCD phase transition, one has /?vac — 42.36k and /vac = (34/3)/51.25 ~ 
0.2211. If K is of order unity, 6vac — 9.37k is probably large enough to render valid the WKB 
result (|4.12|) . This yields an attraction factor Pqcd — 0.49a^/^K^^/^ cos^. In the case of 
the electroweak phase transition, rough estimates give 6vac — (A/4)k where A denotes the 
quartic self-coupling of the Higgs. Its seems therefore probable that ^electroweak < x, so that 
the electroweak transition has only a weak attracting effect on ip. We conclude that phase 
transitions seem to have only a modest effect on ip. It would be at present meaningless to 
refine the calculation of the effect on ip of the electroweak and QCD phase transitions [even 
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the order of the transitions is in doubt, not to mention the precise redshift dependence of 
\{p) and Svac(p)]- In fact, until one has some understanding of the cosmological constant 
problem, it does not make much sense to compute any gravitational effect linked to phase 
transitions. In the following we shall therefore neglect the effect of the phase transitions 
with respect to that of the nine mass thresholds discussed above. 

V. EVOLUTION OF THE DILATON DURING THE MATTER-DOMINATED 

ERA 

The matter content of the universe near the end of the radiation era and during the 
subsequent matter era can be described as the superposition of a relativistic gas ( "radiation" , 
i.e. photons and three neutrinos in the standard picture) and of a non-relativistic one 
("matter"; made of particles of mass mm{<^))- From Eqs. ( |3.5| ) the source terms for the 
cosmological evolution equations ( p.2| ) read p = Pr+Pm, P = Pr+Pm, cr = —am{'-p){Pm-~^Pm) 
with Pr = Pr/S, Pm — 0, and am{^p) = i91n mm{^p)/dip. Either from the definition ( |3.5a] ) 



or from the separate energy balance equations discussed below Eq. (|3.6| ), one deduces that, 
during the expansion, pr oc a~^ while pm oc ■mm{y:>)a~^- Finally, the evolution of ip with 
respect to the p-time p = Ina + const, is given by the equation 

y+[l-X{p,^)]<f' = -[l-3X{p,<f)]aUf), (5.1) 



3 - (/?'2 

with 3X{p,ip) = [1 + Cmm{'^)e^]~^ , C being some constant. In the approximation where 
the radiation era has already attracted ip very near a minimum ipm of mm{<p), we can 
consider that mm{<p) — const, in X{p,ip). Choosing now the origin of p at the equivalence 
between radiation and matter [pr{p = 0) = Pm{p = 0)], we get simply X{p) = 3~^(1 + e^)~^. 
Neglecting (p'"^ in Eq. (|5.1|) and using the harmonic approximation ([4.6|) , we find that (p 
satisfies a linear differential equation which can be rewritten as a hypergeometric equation. 
Denoting x = e^ = a/agquivaience we have 

X{X + l)dlip + (^^X + 2) 9,(^ + ^PmiV -Vm) = 0. (5.2) 

The condition of regularity of ip when x ^ 0, say ip{x = 0) = y^rad (v^rad denoting the value 
of ip at the end of the radiation era, before the transition to the matter era around p = 0), 
selects uniquely the solution of ( |5.2|) to be ipm + (v'rad — Vm) x P[(i, b, c; —x]. Here F[a, b, c; z] 
denotes the usual (Gauss) hypergeometric series. The values of the parameters are 

3 3 

a = - — iu , b = - + iij , c = 2, (5.3) 



with u = 



[fim — I ) -In other words, the attraction factor of the matter era up to the 



1/2 
present time, F^ = {(Pnow - ^m)/{^rs.d - 'Pm), is given by 

F^ = F[a,b,c;~Zo] , (5.4) 



where Zq = e^° = anow/c^equivaience dcuotcs the (Einstein frame) redshift separating us from 
the moment of equivalence between matter and radiation. As Zq is large (see Eq.( |6.5aD 
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below), we can use the asymptotic behavior of the hypergeometric function (together with 
the properties of Euler's F function) to get more exphcit forms for F^- Whatever be the 
sign of Pm — 3/8 (i.e. in the two cases where u is real or pure imaginary) one can write 

with r2 = T{2iuj)/T{2iu + 3/2). Actually, from the estimate ( [4.6b| ) we expect Pm to be 
(much) larger than 3/8 (indeed, if rrim ~ 1 GeV, one would need k to be smaller than 
9.2 X 10^^ to make (3m < 3/8). In that case {u real), one can compute the modulus of the 
complex number r2 in terms of elementary functions to get 



cotanh(27rco') ^ 
7iu{u^ + 1/16) 



fPO 



cos^o ) (5.6) 



with ^0 = ^Po + 2^17 In 2 + Arg(r2). As in the case of attraction by mass thresholds (when 
PaJa ~ 1)5 the attraction factor ( |5.6| ) is proportional to a cosine (when Pm > 3/8) because 
Eq. ( |5.1|) describes a damped oscillation around the minimum ipm of lnmm{<^)- When 
Prn < 3/8, (p slowly rolls down toward ipm without oscillating (overdamped oscillator). [See 
also Ref. |jl3|] in which the transition between radiation domination and matter domination 



was approximated — in the analytical formulas — as being a sharp one] . 



VI. OBSERVABLE CONSEQUENCES OF A COSMOLOGICAL RELAXED 

MASSLESS DILATON. 

Sections 4 and 5 have exhibited several efficient mechanisms for driving the VEV of the 
dilaton toward a value where it decouples from matter. However, none of these mechanisms 
is a perfect attractor. The important question remains of giving quantitative estimates 
of the residual coupling strength of the dilaton at various cosmological epochs and of the 
corresponding observable effects. 

The quantitative estimates of the efficiency of the cosmological attraction of the dilaton 
depend very much on the universality, or lack thereof, of the dilaton couplings. If the dilaton 
coupling functions Bg{^), i?j?($), -Bh($) (the latter representing the class of couplings to 
the fundamental Higgs sector, if it exists as such) are unrelated functions, one expects the 
mass functions ( |2.10| ) to have minima (if any) at different values of ip, say y?^. For instance, 
the lepton and quark masses will involve Bh while hadron masses will all be proportional 
to Bg^^"^ ex-pl—Sn'^b^^ksBs] {B^ = Bsu{3))- In such a non-universal case, the various mass 
thresholds, and phase transitions, will not attract (f to the same value, but will tend to 
reshuffle each time the value of ip. In that case, the only efficient fixing of the value of (p 
would arise during the matter era, ip being attracted toward of minimum of rrimi'p) where 
the label "m" represents the type of matter which dominates the present universe. 

By contrast, one can consider the case where all the dilaton coupling functions coincide, 
i?j($) = B{^). This case of universal coupling of the dilaton to matter has a suggestive 
simplicity. It looks like a natural generalization of the universal e~^* coupling arising at the 
string tree level. In the universal -B($) case, all the mass thresholds, as well as the QCD 
phase transition and the matter era, tend to attract (p to a. common value, some maximum 



12 



(fm of -B ((/)). In the universal case, the cosniological evolution is an extremely efficient way 
of pinning down the value of (p. Moreover, as y? is pinned down to an extremum of B{(p), 
i.e. to a value where dB{(p)/d(p and dmA{(p)/d(p vanish, one can say that the universal 
dilaton coupling case illustrates some "Principle of Least Coupling" in the sense that the 
universe is attracted to dilaton values extremizing the strengths of the interaction. It would 
be worth exploring whether imposing this universality provides a sensible way of selecting 
a preferred class of string models. In the following, we leave open the two possibilities, 
universal/non-universal, in our discussion of the observable consequences of our scenario. 

The earliest observational information we have about cosmology concerns the primordial 
abundance of the light elements (mainly Helium 4, with traces of Deuterium, Helium 3 and 
Lithium 7). Let us discuss the production of Helium 4 as an example. In the standard 
scenario of homogeneous primordial nucleosynthesis, the abundance of Helium is mainly 
determined by the neutron/proton ratio at the temperature where the rate of interconversion 
n <->■ p due to weak interactions becomes slower than the cosmological expansion rate (freeze- 



out) (see Ref. |jT9|). Neglecting the small additional effect of free neutron decay, one can 
write an approximate analytical formula for the primordial Helium abundance (by weight), 
Y, of the form, Y = 2/(exp(aX) + 1) where a is a pure number of order unity and where X 
denotes the following dimensionless combination of coupling constants and masses 

X . /Ai + sAr' (!^i^^) r^^S^y" {-ir-y'" . (e.D 

V mw / V mw / V 10.75/ 

Here g2 denote the SU{2) coupling constant, qa — 1-26 the axial/vector coupling of the 
nucleon, and g^ the effective number of relativistic degrees of freedom at freeze-out (retained 
here to allow easy comparisons between the effect of a change in g^, — e.g. an additional light 
neutrino — and the effects of changing, e.g., Newton's constant G = mp^^^^^^, or Fermi's one 
Gp = 5'2/8'^w/)- A remarkable fact about the combination X is that it is numerically of order 
unity thanks to a dehcate compensation between large {mpia_nck/mwY^'^ — (1-52 x 10"'^'^)"'^/^ 
and, small [(m„ — nip) /my/ c^ L61 x 10^^] factors. This fact prevents us from proposing an 
educated guess of the quantitative dependence of X on the bare dilaton coupling constants 
Bf{i~p)i Bg{(f),. . . . Even the sign oi d In X / d\n B~^ (when Bi{ip) = B{(f)) is unclear. On the 
other hand one can estimate that dY/dlnX ~ —0.44 both from the rough analytical formula 
for Y{X) and from the numerical computations of the dependence of Y on the neutron half- 
life or on g^:. We can therefore write the value of the Helium abundance predicted by a 
scenario modified by the presence of a dilaton as 

Fdii(^) = yGR(^) _ 0.22^^«:((/.,,d - Vmf , (6-2) 

where we have reestablished the slight dependence of Y upon the baryon to photon ratio, rj. 
In the standard, general relativistic scenario the dependence of the GR-predicted abundances 



^Note that a primordial (inflationary type) phase transition — as well as the electroweak one, 
if the Higgs sector is fundamental — could instead attract ip to a minimum of B{ip) through a 
transient vacuum energy oc B{ip). 
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on T] is crucially used, together with the observed values of the light-element abundances, to 
set upper bounds on t], and thereby upper bounds of the ratio of the present total baryon 
mass density to the closure density, fih- The standard conclusion being that baryons fail 
to close the universe by at least a factor five, flh < 0.2 |19|. Eq. ( |6.2|) [to be completed by 



the corresponding dilaton-modified predictions for the other light elements] suggest that a 
dilatonic universe could naturally accomodate Qb = 1 ii the value yjrad of ip at freeze-out (i.e. 
just after the electron mass threshold) differs by a small (but not too small) amount from 
the minimum ip^ [For instance, in the case of the Helium abundance, the dilaton correction 
term on the right-hand side of Eq. ( |b.2| ) should be approximately —0.03, and dlnX/dlnB~^ 
should be positive] . It would be interesting to reexamine in full numerical detail primordial 
nucleosynthesis within the type of dilaton scenario considered here to assess whether it 
could naturally reconcile Qb = ^ with the observed abundances of light elements. Let us 
only note here that the rather modest attraction toward (pm which is probably needed in 
such a scenario seems more natural in the non- universal case. Indeed, in the universal B{(p) 
case, all the nine mass thresholds compound their effect to drive (p very near some universal 
minimum (pm- More precisely, yjrad — 'Pm = Fr x (ipi^ — iprn) where ipin is the "initial" value 
of ip (meaning in this work, before the electroweak phase transition) and where the total 
attracting power of the radiation era is given by 

Fr{K) = I n m^iPAfj)] X I n piprin \ ■ (6.3) 

{a=i ) [4=2,3 J 

The values of jSa and /™ to be used in the attraction factors of each of the nine mass 
thresholds have been given above [remember that the +(— ) sign corresponds to fermions 
(bosons)]. The second factor in Eq. ( |0|) corresponds to the effect of the two known phase 
transitions electroweak (2) and QCD (3). In view of the uncertainty in the calculation of the 
effect of phase transitions, and anyway of their expected modest contribution (see above), 
we shall neglect the attraction power of these phase transitions in the following. The small 
but non-zero value of (pra.d — 'Pm = Fr{K)A(p [with Aip = ipi^ — ipm] implies that all the gauge 
coupling constants squared, g"^ oc i?~^ (</)), differed, at the end of radiation era, from their 
present values g^ by a fractional amount 



^'^ , ^° ^ ^«:(^rad - Vn.f = i«:(i^,(«:) A^)^ (6.4) 

[where we used the fact that (p^ — (pm ^ 'PiaA — 'Pm because of the matter era attraction]. 
As one a priori expects Aip = (p^^ — ip„i to be of order unity, the function )^kF'^{h), which is 
plotted in Fig. 2, illustrates the remarkable efficiency (in the universal case) of the radiation 
era in pinning down the values of the physical coupling constants. 

During the subsequent matter era, p) is (in the universal case) further driven toward y^m 
by the factor Fm{n, Zq), Eqs. (|5.4|) - (|5l^) . The numerical value of the matter-era attraction 



factor Fm is proportional to Zq where Zq = e'^^ denotes the redshift separating us from the 
epoch of equivalence between matter and radiation. In the approximation mm{<p) — const, 
introduced at the beginning of Sec. 0, this redshift is given by ||T3| 



^matter 

^0 = ^^^T^ - 13350 fir5 , (6.5a) 



^rad 

Po 
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where 

^75 = 87rGpr"'V[3(75 km s^^Mpc"^)^] = pr"'7l-0568 x lO'^^g cm-=^ . (6.5b) 

Under the hypothesis of a spatially flat universe {K = 0), generally assumed in 
this paper, Q-r^ is linked to the present value of Hubble's "constant", Hq, by ^75 = 
{Hq/75 km s~^Mpc~^)^. In that case, the observational limits 50 < Hq/1 km s~^Mpc^^ < 
100 imply 0.44 < ^^75 < 1.78. On the other hand, if one assumes that the universe is spa- 
tially hyperbolic {K = — 1), one must modify the coefficients of the evolution equation ( p.3[ ) 
for ip by retaining the i^-dependent terms. However, it was shown in |T^ that as long as 



f275 > 0.05 this modification of Eq. (|3.3|) has a small effect, and that the matter-era attrac- 



tion factor oi K = —1 universes is well approximated by the K = formula ( ^.6|) , with Zq 
given by Eqs. ( |6.5aD , ( |6.5b| ). The main difference is that now ^75 is not related to Hq, and 
can be smaller than 0.44. In fact, present observational data are compatible with ^75 ~ 0.1. 
Finally, the scenarios considered here predict that the present value of ip, say ipo, differs 
from the minimum (pm by (po — (pm = Ft{K, Zo)Aip where Aip = y^in — ipm and where the 
total attraction factor is 

FtiK, Zo) = FriK) Frr^iK, Zo) . (6.6) 

There are three kinds of presently observable consequences of having ipo near, but different 
from, (fj^: (i) violations of the (weak) equivalence principle; (ii) modifications of relativistic 
gravity, and; (iii) slow changes of the coupling constants of physics, notably the fine-structure 
constant a and Newton's constant G. 

To discuss the modifications of the gravitational sector, we can make use of the re- 
sults of Ref. 1^ on the relativistic gravitational interaction of condensed bodies in generic 



metrically-coupled tensor-scalar theories. Indeed, the action describing the classical inter- 
action of massive particles of various species under the exchange of the g^u and ip fields is 
given by Sg^^ + Sm[g, 'P, x] where Sg^^, is given by ( |2.7b| ) and Sm by Eq. (|3.4|). This action is 



identical to the one studied in Sec. 6 of pOf. We conclude that, at the Newtonian approx- 
imation, the interaction potential between particle A and particle B is —GABfT^AfT^B/^AB 
where tab = |xa — x^l and 

GAB = G{l + a^2'^aP) . (6.7) 

Here G is the bare gravitational coupling constant entering the action ( p.7b| ), and a}^ is 
the present strength of the coupling of the dilaton to A-tjpe particles, i.e. the value of ( p.6|) 
taken at the cosmologically determined VEV (po. [In diagrammatic language, the two terms 
on the right-hand side of Eq. ( |6.7| ) are, respectively, the one-graviton exchange contribution 
(G) and the one-dilaton exchange one {Ga]^ a^ )]. Two test masses, made respectively of 
A- and 5-type particles, will fall in the gravitational field generated by an external mass 
niE with accelerations ua and qb differing by 

Ax f (0) (0)x (0) 

a J AB ttA + ttB l + i(a7 + "B)"E 

All precision tests of the gravitational interaction used macroscopic bodies made of (neutral) 
atoms. Let the labels A, B,. . . denote some atoms. In the approximation where one neglects 
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'fnuI'mN, "^dl'^Ni 'raf^jm^^ a and a„cak, the mass of an atom is a pure (dilaton-independent) 
number times a QCD-determined mass scale, say u^{'~p). In this approximation aA(v^) = 
d In mA/dip is independent of the type of atom considered and is equal to a^.^^) = 5 In u^^jd^. 
The dilaton dependence of Uj, is determined by Eq. ( p.9|) . Choosing uj, so that its present 
value uj,{i{)Q) is numerically equal to the atomic mass unit, u = 931.49432 MeV, we see from 
Eqs. (H^) that 

aA(v5o) ^ Oisiipo) = j3-i{ipQ - ipra) = PsFti^, Zo)Alp , (6.9) 

with p3 = Kdlnus/dlnB-^ ~ 40.82 /t. 

In this approximation, the dilaton mimics a usual Jordan-Fierz (-Brans-Dicke) field, i.e. 
a scalar field coupled exactly to T/^. The main observational consequences of the body- 
independent coupling ( p.9|) are modifications of post-Newtonian relativistic effects, 0{v^ /(?) 
beyond the Newtonian 1/ R interaction (weak gravitational field case )0. The latter are 
measured by the two Eddington parameters 7Edd — 1 and /^Edd — 1 (which vanish in general 
relativity). From |2^ we see that in the approximation (|6.9| ) 



1 - 7Edd = 2-^ ^ 2(/33)'(Fi(«:, Z^)A^f , (6.10) 

/^Edd - 1 = \t^^ ^ \m\F,{K, Zo)Avf . (6.11) 

2(l + a|) 2 

Note also that the value of Newton's gravitational constant (in Einstein units) is G^ = 
G{l + al). 

Much more sensitive tests of the existence of dilaton couplings are obtained by looking at 
violations of the weak equivalence principle, i.e. at the body-dependence of aA{^o) beyond 
the QCD approximation (|6.9| ). To do this, we shall retain the leading 771^/711^, md/rrij^, 
me/rriN and a corrections to the mass of an atom. First, the mass of the nucleons have 
the form, rup = niNs + burriu + bdrud + CpO, m„ = 771^3 + hdirtu + buTUd + CnOt-, where 
mAr3(^ M3) is the pure QCD approximation to the nucleon mass, and where &„, 6^, Cp/u^ 
and Cnjuj, are pure numbers (in the approximation of negligible strange-quark content one 
has hu = {p\uu\p)/27nN, bd = {p\dd\p)/2mj^) []22[. Second, the mass of an atom can be 
approximately decomposed as 

m(Atom) = Znip + Nm^ + Zm^ + Eg^"^'^"^ + ^""^1^"^ , 

where Z is the atomic number and A^ the number of neutrons, and where E!^^'^^'^^^ denotes 
the strong-interaction contribution to the binding energy of the nucleus, and i^p^^'^"'^ the 
Coulomb interaction energy of the nucleus. 

In terms of the baryon number B = N + Z, the neutron excess D = N — Z, and the 
Coulomb energy term E = Z[Z — 1)/{N + 2')^/^, the mass of an atom can be written as. 



^In view of the positiveness of /33, the recent results of [21| show that the deviations from general 
relativity are further quenched in the strong- gravitational-field case of binary neutron star systems. 
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m(Atom) = M3M3 + a'B + 5'D + a^ausE , (6.12) 

where M3 is a pure number (= B+ strong-interaction binding contribution) and where we 
have defined 

111 1111 

a' = a + -Cna + -CpO + -m^ , 6' = -6 + -CnQ - -CpO - -rUe , 

with the usual definitions for a = \{mu + md){hu + hd), 5 = {rrid — mu){hu — hd). [Note the 
factor 1/2 in the first term of the definition of 5']. Finally, by differentiating the logarithm 
of ( |6.12| ) we get a more precise expression than ( |6.9| ) for the dilaton coupling strength 

38 f^. ^^/m ^„,^/^^ 



a.Cpo) - a,W.) + a^ (m j, + «^ (m j, + «3g^ [j-^)^ , (6.13) 

where we have introduced a = a'/u^, 6 = 5' /u^, and approximated M3 ^ M = m(Atom)/M3 
in the corrections terms.Q Finally, from Eq. ( |6.8|) we get an equivalence-principle violation 
of the form 



— ) ={KFt{K,Z^)^^f 

a / AB 



^-^Is) + ^"^(f)+^-^(s 



AB 



(6.14) 



where {AX)ab = Xa — Xb and where Cb = Xuid^ /d\nB^^ , Cd = Xu^dS/dlnB^^, Ce = 
XugXadsa, A„3 = dlnus/dlnB'^ and A^ = dlna/dlnB^^. Numerically, our usual estimate 
( [2. 91 ) gives Aug ~ 40.82, and the idea of unification of gauge couplings at the string scale 
gives Aq ~ 1. [E.g. in the simplest S'?7(5)-type GUT the value of the fine structure constant 
at the QCD-confining energy scale M3 — such that astrong(^3) — 1 — is given by 

a{us)-' = (22/7)aG^T - (10/21)astrong(%)'' ^ (22/7)aGUT ^ B{^)] . 

We have also a^a = 0.717 MeV/M3 = 0.770 x 10~^ from the fit of atomic masses to the 
Bethe-Weizsacker formula. We can therefore estimate the coefficient of the nuclear Coulomb 
energy term in Eq. (|6.14 ) to be Ce — 3.14 x 10~^. As for the other two coefficients, Cb 



and Cn, it is much less clear how to estimate them. From the experiment-derived values of 
(T = 35 ± 5 MeV and 6 = 2.05 ± 0.30 MeV, plus the theoretical estimates Cpa = 0.63 MeV, 
C„a = -0.13 MeV P|, one can compute a = 3.8 x 10"^ and J = 4.2 x 10"^ From the 



point of view of their dilaton dependence a and 6 are the sum of four terms proportional 
to niu/us, rrid/us, rUe/u^ and a. It is impossible at present to rehably guess the dilaton- 
dependence of the mass ratios mquark/^^hadron and me/mhadron- The numbers we would get 
for da/dlnB^^ would be very different were we to assume our usual exponential link to 
the string scale, or some other assumption. It seems however reasonable to estimate that 



^Note that the assumption of a universal B{ip) is crucial to ensure that all the terms in Eq.( |6.13 ) 



have in common a very small factor ipo — ipm- If, e.g., the mass of leptons (and/or a) depended 
on a different function of 99 than the mass of hadrons, Eq.( |6.13| ) would, at best, predict that the 
equivalence principle is violated at the (unacceptable) level ~ {0{a) + 0(w.iepton/"iN))^, in the 
favourable case where the universe is assumed to be dominated by hadronic matter. 
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the order of magnitude of da /din B^^ and d5/dlnB~^ will be at most that given by the 
exponential assumption ( |2.11| ), and at least that obtained by differentiating only the fine- 
structure constant contributions to a and 6. This yields corresponding rough upper and 
lower bounds for the coefficients of the B and D contributions: 1.1 x 10~^ ^ \Cb\ ~ 7.5, 
1.7 X 10^^ ^ \Cd\ ~ 8.2 X 10^^. If these upper bounds are correct, one can check that 
the last term in Eq. (|6.14|) will be numerically dominant for pairs {A, B) having a large 



difference in atomic number [Indeed, E/M is roughly proportional to Z'^^^]. The largest 
effect would arise in comparing Uranium (for which E/M ~ 5.7) with Hydrogen (any light 
element would do nearly as well). For such a pair, Eq. ( |S.14| ) yields 

(Aa/a)^ax = O.lSinFtiK, Zo)Aipf . (6.15) 



The right-hand side of Eq. ( |6.15| ) is plotted in Fig. 3 as a function of k (assuming Aip = 1, 



and (^75 = 1). As one a priori expects k to be of order unity. Fig. 3 shows that, within 
the scenario considered here (including universal dilaton couplings), the present tests of the 
equivalence principle (at the 10~^^ — 10~^^ level) do not put any significant constraints on 
the existence of a massless dilaton. 

The situation is even worse if we consider tests of post-Newtonian gravity. Indeed, from 
Eqs. (|6.10| ) and (|6.15|) we have the link 

""^ 5.4 X 10-^(1 - 7Edd) (6.16) 



showing that the present and planned levels of testing of post-Newtonian gravity, i.e. 10 ^ 
and 10~^ at best for 7Edd, correspond, respectively, to equivalence-principle tests at the levels 
5 X 10"^ and 5 x 10"^^. The other link /?Edd - 1 = iPsi^ - 7Edd) ^ 10.2fi;(l - 7Edd) shows 
that /^Edd tests do not fare essentially better. 

The last observational consequence of our scenario to discuss is the residual present 
variation of the coupling constants of physics. From dlna/dlnB^^ ~ 1, with InB^^ = 
const. + ^K,{(p{t) — (fmY and a present time dependence of ip given by Eq. (|5.6| ) [in which 
the leading term is cos{ujp + const.)], we deduce that 



aH 



a 



3 

cjtant^o H — 



{Ft{K,Zo)A^f . (6.17) 



Using u = [|(/5m — |)]^^^ with Prn ~ 40.8k (if the particles dominating the universe have 
a mass not very different from the GeV scale), we have the approximate link {a/aH)o ~ 
— 43k~^/^ tan6'o(Aa/a)niax- This link shows again that equivalence principle tests are the 
most sensitive way of searching for possible dilaton couplings [In the foreseeable future, 
ultrastable cold-atom clocks might probe the level a/a ~ lO^^^yr^^ ~ IO^^Hq which corre- 
sponds to Aa/a ~ 10^^]. 

Finally, the time variation of the gravitational coupling constant (measured in Einstein 
units) is 
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{PsFt{n,Zo)A^f . (6.18) 



In actual G experiments one is comparing an orbital frequency n (e.g. let us consider that 
of a planet around the Sun) to an atomic frequency z/. Taking into account the adiabatic 



invariants of the orbital motion (angular momentum and eccentricity) and assuming an 
atomic clock based on the Bohr frequency oc mf.a'^, the directly measured quantity will be 

(6.19) 
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where rris, nip denote the masses of the Sun and of the planet. [Contrary to Eq. ( |6.18D 
Eq. ( |6.19| ) is valid in any system of units] . Eq. (|6.19| ) gives finally^] 
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(Fi(«:,Zo)Ay.)2 4/31 + 5/53 -/3e-2«: , (6.20) 



in which the term coming from Eq. ( p.l8| ) dominates. In spite of the large factor {j3^/ nf — 
(40.8)^, the scaling of the prediction ( |6.20| ) with the Hubble rate Hq makes it pale in com- 
parison with equivalence principle tests. [On the other hand, this large factor renders G 
experiments competitive with a ones]. 



VII. CONCLUSIONS 

Einstein's starting point in constructing general relativity was the interpretation of the 
universality of free fall in terms of a universal coupling of matter to a common metric tensor 
g^j,y. It has since been felt that such a universal metric coupling was the only theoretically 
natural way of explaining how the long-range fields participating in gravity]^ could satisfy 
the high-precision tests of the equivalence principle (now reaching the 10~^^ level). 

The present work suggests that a universal multiplicative coupling of a long-range scalar 
field $ to all the other fields, £tot = B{^)CQ{g^ui ^^ ^m, ^5 • • •) with -B($) admitting a local 
maximum, though a priori entailing strong violations of the equivalence principle, provides 
another theoretically natural way of explaining why no violations have been seen at the 
10^^^ level. It maybe worthwhile to summarize in qualitative termsj^ the basic reasons why 
a massless dilaton is rendered nearly invisible during the cosmological evolution: (i) Each 
time, during the radiation era, the universe passes through a temperature T ~ tua the 
A-type particles and antiparticles become nonrelativistic before annihilating themselves and 
disappearing from the cosmic soup; this provides a source term for the dilaton proportional 
to the (y9-gradient of ?7i^ (</)), which attracts ip toward a minimum ip^ of ?7i^ (</)); Eq. (|4.6b[) 



and Fig.l suggest that each such attraction is moderately efficient, leaving Lp nearer to y?^ by 



^The link between ipQ and ii/n is more involved if n is the orbital frequency of a binary neutron star 



system: see [23|, which must be completed by taking into account the changes in the rest-masses 



of the stars, and the non-perturbative gravitational self energy effects pi[ . 

^In Einstein's theory gravity is mediated by only one, spin 2, field; but in metrically-coupled 
tensor-scalar theories, gravity is mediated both by a spin 2 and a spin field. In the latter case, the 
universal metric coupled to matter is a combination of the two pure-spin fields: 5!^°'^ = A^{ip)g* . 

^*^Simplified quantitative estimates are provided below. 
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a factor ~ 1/3. (ii) Under the assumption of universality of the dilaton couphng functions 
B{(f), the minima of all the mass functions mAi^p) will coincide and the ~ 9 mass thresholds 
of the radiation era will compound their effects to attract very efficiently (f toward some 
common minimum ip^- (hi) In the subsequent matter era, (p will be continuously attracted 
toward a minimum of the mass function rrimi'p) corresponding to the (nonrelativistic) matter 
dominating the universe. [ Under the same universality condition this minimum will be again 
(pm-] The attraction factor due to the matter era is inversely proportional to the 3/4th power 
of the redshift Zq ~ 1.3 x 10^ separating us from the end of the radiation-dominated era. (iv) 
As a consequence of the very efficient total attraction toward ipm, the present strength of the 
coupling of the dilaton to any type of matter aAi'p) , being proportional to the (/^-gradient 
of m^ ((/)), is very small. The present deviations from general relativity in the interaction 
between two masses, uia and tub, are proportional to the product aAC^B and are therefore 
extremely small, (v) The equivalence principle tests are very sensitive, but they probe only 
differences {a a — otB)ac which, because of the known universal features of QCD-generated 
masses, contain as supplementary small parameters either the ratio of the quark masses to 
the nucleon mass, or the fine-structure constant. 

From a theoretical point of view, our work suggests a criterion for selecting a preferred 
class of string models: namely those where string-loop effects preserve the universal multi- 
plicative coupling present at tree-level, with a dilaton-dependent function admitting a local 
maximum^. It will take, however, an improvement in our current understanding of super- 
symmetry breaking in string theory to see whether the universality required by the Least 
Coupling Principle is a viable option, providing a reasonable selection criterion for SUSY 
breaking mechanisms. It is to be noted that in this paper we had always in mind the cou- 
pling of the (four dimensional) dilaton which is such an intimate partner of the graviton 
that it seems reasonable to assume that it remains massless in the low-energy worlcQ. How- 
ever, the cosmological attractor mechanism described here could also apply to the other 
gauge-neutral scalar fields (moduli) present in string theory. Because of threshold effects, 
the gauge coupling function Bp acquires a non trivial dependence on the moduli fields p4| . 



Therefore, we have here a possible mechanism for fixing the moduli to values where they 
decouple from the other fields. 

From an experimental point of view, our results provides a new incentive to improving 
the precision of equivalence principle tests (universality of free fall, constancy of the con- 
stants,. . . ). Fig. 3 suggests, when assuming that the curvature of B{lp) near its maximum 
is of order unity — say 0.1 < k < 10 — , to look for a present level of violation of the uni- 
versality of free fall somewhere between 10~^^ and 10^^^. Actually, one should not consider 



^^The existence of a local maximum in -B($) is necessary. For instance, in the case of the tree- 
level coupling function exp(— 2$), $ would continuously roll toward — oo during the cosmological 
expansion, and worse, by Eq. (^^), $ would cause deviations from general relativity, including 
violations of the equivalence principle, of order unity or more. 

^^It is enticing to assume that the presently obscure mechanism ensuring the vanishing of the 
cosmological constant allows both 'gravitational' fields to remain long-ranged in the low-energy 
world. 
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the results plotted in Fig. 3 too seriously. On the one hand, even within the precise assump- 
tions made in the text, the predicted maximal value of Aa/a contains an unknown factor 
~ ^^75 {(pin — (pm)'^ which could be ^ 10. On the other hand, our assumption ( p.llD with 



/iA = 1 has entailed a specific phasing of the various oscillations undergone by if during the 
radiation era, i.e. specific choices of where on the curves of Fig. 1 {ip{p) — iprn)/{'f{—oo)—(prn) 
ends up being when p — > +oo. It is possible that the assumption (p.ll|) has overestimated 
the combined attraction power of the radiation era mass thresholds. A different estimate 
is obtained by multiplying the WKB approximations ( [4.10a| ) (all valid as soon as k ^ 1), 
assuming that all the oscillation angles 6^ are randomly distributed on the circle. Under 
the latter assumption it makes sense to compute a rms value of the radiation era attraction 
factor ((cos 0^)rms = 1/ v2)- Neglecting as above the effect of the phase transitions one finds 
[Fr{K)]rms = 1-87 X 10~^ X K"^/^ for the radiation era attraction factor, and from Eq. ( p.6| ) 
with cu 3> 1, {cos9o)rms = 1/v^ and ^75 = 1, [Fm{K)]rms = 1.47 X 10~^ X K~^/^ for the matter 
era attraction factor. This leads to a total attraction factor [Ft{K,)]rms = 2.75 x 10~^ x k~^ 
and to the following analytical estimate of the rms value of the maximum value of the 
equivalence principle violation 

{Aa/a)Z: = 1-36 x 10-i««:-^(A^)2 . (7.1) 

The comparison of Fig. 3 with Eq. (|7. 1|) (valid if k ^ 1 and the angles 9-^ are randomly dis- 
tributed) indicates that the phasing of the radiation era oscillations tends to be destructive. 
It is possible that alternative assumptions, different from ( p.ll| ), yields values of (Aa/a) 
nearer to the rms analytical estimate ( [r.l|) . This would have the consequence that presently 
planned satellite tests of the equivalence principle [^ which aim at the level Aa/a ~ 10~^^, 



would probe a larger domain of values of k, Aip and ^75. 

In conclusion, high-precision tests of the equivalence principle can be viewed as windows 
on string-scale physics. Not only could they discover the dilaton, but, by fitting observed 
data to the expected composition dependence ( |6.14| ) of the equivalence principle violation, 
they could give access to the ratios Cb/Ce, Cd/Ce which are delicate probes of some of the 
presently most obscure aspects of particle physics: Higgs sector and unification of coupling 
constants. 
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FIGURES 

FIG. 1. The factor by which 99 is attracted (when the early universe cools down through 
T r^ ttia) toward a minimum ipm of the function rriAi^) is plotted as a function of b^ = PaJa- 
The solid (dashed) line corresponds to A being a fermion (boson). 

FIG. 2. The solid line represents logio[(fl'rad~5o)/fl'o] ^^ ^ function of log^g '^j i-^- the fractional 
deviation (left over at the end of the radiation era) of the gauge coupling constants g'^^^ oc i?~^((/^rad) 
from their present values ^g, versus the curvature k of the function lnB~^{ip) near its minimum. 
The dashed line represents an analytical estimate (when k ^ 1) of that deviation, obtained by 
assuming that the phases 9 of the WKB results Eq. ( |4.10a|) are randomly distributed. 

FIG. 3. The solid line represents logio(^'^/'2)max as a function of log^^Q k, i.e. the expected 
present level of violation of the equivalence principle (when comparing Uranium with a light el- 
ement) as a function of the curvature k of the (string-loop induced) function lnB~^{ip) near a 
minimum ipm- The dashed line represents an analytical estimate (when k ^ 1) of that violation 



obtained by assuming random phases 9 in Eqs. (4.10a) and (|5T 
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